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r^ . Abstract. Local Scale Invariance (LSI) is a theory for anisotropic critical phenomena 

r^ I designed in the spirit of conformal invariance. For a given representation of its 

generators it makes non-trivial predictions about the form of universal scaling 
"jrt . functions. In the past decade several representations have been identified and the 

"^ ' corresponding predictions were confirmed for various anisotropic critical systems. Such 

tests are usually based on a comparison of two-point quantities such as autocorrelation 

and response functions. The present work highlights a potential problem of the theory 

in the sense that it may predict any type of two-point function. More specifically, it is 

i-rt ■ argued that for a given two-point correlator it is possible to construct a representation 

SIh ' of the generators which exactly reproduces this particular correlator. This observation 

^ ■ calls for a critical examination of the predictive content of the theory. 
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Q ' 1. Introduction 

'k>( I Local scale invariance (LSI) stands for a theory developed by M. Henkel and 

H I collaborators which generalizes global scale invariance of anisotropic critical systems 

and ageing phenomena to a local space-time-dependent symmetry [1]. It is inspired by 
the success of conformal invariance applied to two-dimensional equilibrium systems [2,3] 
and uses a similar terminology. 

Local scale transformations for anisotropic systems are generated by an infinite set 
of generators X_i, Xq, Xi, X2, . . . and Y_i/z, Yli/^+i, Y^i/z+2, ■ ■ ■■ These generators obey 
the commutation relations 



[Xn, Xm] = {n- m)Xn+m (l) 

where z is the dynamical exponent which quantifies the degree of anisotropy 
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Remarkably, the theory predicts the specific form of universal scaling functions 
appearing in response or correlation functions while it makes no prediction about 
the values of critical exponents. In recent years the theory has been extended and 
successfully applied to a large variety of models [4-21]. This success raised the hope 
that LSI could be a generic symmetry of scale-invariant anisotropic systems. However, 
some authors reported results which seem to be incompatible with the predictions of 
LSI [22-29] which released a debate concerning the applicability of the theory. 

In a given model local scale invariance can be established by choosing a suitable 
representation of the algebra dl])-® in such a way that the system under consideration 
is invariant under transformations generated by X„ and Y^. In some cases such a 
representation can be derived exactly from an underlying partial differential equation, 
whereas in other cases LSI is just assumed as a hypothetical symmetry, leading to 
certain predictions which can be tested by numerical simulations. Since LSI is a model- 
independent theory these predictions depend exclusively on the chosen representation 
of its generators. 

The simplest representation can be derived from the Schrodinger equation which 
describes diffusing particles with z = 2 [1]. This representation was generalized to the 
case z =/^ 2 and successfully testet for z = 4 [21]. Recently, Baumann and Henkel [30] 
found two representations of the LSI algebra ([ll)-(l2]) for arbitrary z which involve non- 
local fractional derivatives [31] (see Appendix B). The present work generalizes this 
concept even further by considering representations with arbitrary non-local operators. 
It is shown that the resulting representations of the LSI algebra are so general that 
in principle any two-point correlation function can be reproduced by the theory. This 
has important consequences regarding the predictive power of the theory, as will be 
discussed at the end of this paper. 

2. Geometrical interpretation of the generators 

To gain some intuition how the generators Xn and Ym work let us first recall the 
geometrical interpretation of the generated transformations. For simplicity, we will 
restrict to the l-|-l-dimensional case. The generalization to higher dimensions is not 
difficult and requires to replace Ym by a vector operator Ym ■, see Refs. [1,32]. 

The simplest representation of the generators X„ and Ym-, which describes the 
geometrical content of local scale transformations, is given by 

Xn = -f'+^dt — ^erdr, (3) 

z 

n_l/, = - t^dr , (4) 

where we have adopted the convention of non-integral indices m = k — 1/z with k eN. 
One can easily verify that these operators satisfy the commutation relations ([1]) and ([2]). 
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Figure 1. Distortion of a square lattice in 1+1 dimensions through local scale 
transformations generated by various LSI generators, here shown for the case z = 2. 
Qualitatively similar deformations are observed for any z > 1. 



Moreover, one can see that they carry the physical dimensions 

[Xn] = [time]" (5) 

[n_i/J = [timel'^Pength]-^ (6) 

As shown in Fig. [1], each of these generators corresponds to a well-defined geometrical 

transformation in space-time. For example, the generator X_i = —dt generates 
translations in time while YLi/^ = —dr generates translations in space: 

exp(rX_i)/(t,r) =/(t-r,r), (7) 

exp(sr_i/,)/(t,r) = /(t,r-s). (8) 

Likewise Xq = —tdt — -rdr generates anisotropic dilatations 



exp(AXo)/(t,r) = /(t/6,r/6V- 



(9) 



by the factor b = e^ and thus it can be identified as the generator of global scale trans- 
formations. Here the amount of anisotropy is controlled by the dynamical exponent z. 

The lowest non-trivial operators, which mix space and time, are Xi and l^i-i/^. As 
shown in Fig. [T|, the operator Yi-i/^ generates a Galilei transformation 



exp(cYi_i/^)/(t,r) = /(t, r - ct) 



(10) 



which may be interpreted as a global shear transformation in space-time. Similarly, Xi 
generates a dilatation with an elongation factor proportional to the actual time. By 
combining the action of all generators it is possible to generate the full group of local 
scale transformations. Note that these transformations obey causality, i.e., in Fig. [1] 
horiztonal lines will always be mapped onto horizontal lines. 
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Figure 2. Action of tlie generator Y_i/2 using the example of a simple diffusion 
equation dt(j>{r,t) ~ D'S/'^(f){r,t). In this case the two-point response function is an 

exp(— ^^) which is shown in the left 
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ordinary Gaussian distribution G{r, t) 
panel. As can be seen, this function is not invariant under infinitesimal transformations 
generated by Y1/2 — —td^ since geometrical shear leads to a skewed function of the 
form G{r, i)(l-t-^). Therefore, the generator Y1/2 = —tdr is not a symmetry operator, 
rather it has to be extended by an additional term that compensates this tilt, restoring 
the original form of the function. For the diffusion equation this term takes the simple 
form — A^r with the so-called mass M — 1^. As shown in the right panel, a subsequent 
application of this term restores the original non-tilted function up to order O(e^). This 
demonstrates the mechanism which makes the diffusion equation invariant under the 
action of the generator Y1/2 = —tdr — Mr. 



3. Local scale invariance of quasi-primary fields 



As a next step, the geometric generators ([3])-(j4]) have to be extended in order to describe 
the symmetry properties of a scale-free critical phenomenon. Within the framework 
of LSI it is assumed that the physical properties of such a system can be expressed 
in terms of so-called quasi-primary fields^ here denoted as (f){t,r), which transform 
covariantly under the action of the generators. As we will see, this requires to extend 
the representation of the generators by additional terms. 

For simplicity let us assume that the system under consideration is translational 
invariant in space and time, for example a critical kinetic Ising model in its stationary 
equilibrium state. This means that all quasi-primary fields are translational invariant 
as well, hence the two operators X_i = —dt and Y_i/z = —dr are already symmetry 
operators and need no extension. The situation is different for global dilatations 
generated by Xq. Here a quasi-primary field changes its amplitude according to the 
scaling law 

0(t/6,r/6i/") = 6"/"0(t,r), (11) 

where the exponent x is the so-called scaling index associated with the field (f){t,r). 

I Presently it is not yet fully clear how quasi-primary fields can be characterized in anisotropic systems. 
Usually it is believed that order parameter fields are quasi-primary while their derivatives are not. For 
a discussion see Ref. [15]. 
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Consequently, the bare generator Xq defined in (^ is no longer a symmetry operator, 
rather it has to be extended by a suitable term that compensates the change in the 
field amplitude. Obviously the required term is just a constant, leading to the standard 
representation 

Xo = -tdt - -rdr --. (12) 

z z 

Next, let us consider the operator Fi-i/^ which generates global shear transformations. 
It is important to note that global shear itself is generally no symmetry transformation 
because it distorts the fields in a non-trivial way. This is demonstrated in Fig. [2] using 
the example of the diffusion equation where z = 2. As can be seen, global shear distorts 
the response function, leading to a skewed profile in space-time. Therefore, the generator 
Yi^ijz has to be extended by a suitable additional term that compensates this distortion. 
For the diffusion equation this term takes the particularly simple form —Air, where M. 
is the so-called mass parameter. However, in general the required compensation terms 
may be much more complicated. 

To be as general as possible, we therefore assume that Yi-i/^ is extended by an 
arbitrary linear operator Bi^ijz'- 

Yi-i/z = -tdr - Bi^i/z . (13) 

Likewise we assume that all remaining generators, which generate combinations of shear 
and dilatations, are extended by certain linear operators as well: 

Xn = - r+^9i ^tVS, - -(n + l)t" -An n = 1,2,3,... (14) 

Yk-i/z = - t'dr - Bk-i/z k = 1,2,3,... (15) 

The linear operators An and Bm, by which the generators are extended, are of course not 
independent, rather they are constrained by the commutation relations. For example, 
the operator A2 cannot be chosen freely, instead it is constrained by the commutation 
relations 

[X2,Xo]=2X2 =^ [A2,Xo] = 2A2 (16) 

[[X2,X_i],X_i] = 6Xo =^ [[A,X_i],X_i]=0. (17) 



4. Iterative construction of the generators 

The first point of this work is to show that any space-time representation of the LSI 
commutation relations is fully determined by the generator X2 or, equivalently, by the 
operator A2\3 

The construction starts by choosing the operator A2 in such a way that it obeys the 
constraints ( TT6l) . Once A2 is specified, all other operators can be constructed iteratively 

§ If one is only interested in the subalgebra {X_i, Xq, Xi, y_i/z, ^1-1/2} it even suffices to specify Ai- 
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as follows. The lowest generators X_i, Xq and Vli/^ are always given by their standard 
representation 

X_i := -dt (18) 

Xo := -tdt--rdr-- (19) 

z z 

Y-xiz:= =-dr. (20) 

Moreover, the generator Xi can be computed by setting 

Xi:=i[X2,X_i] (21) 

meaning that Ai = ^[dt,A2]- Likewise, all other generators can be constructed 
recursively by 

^n:=^^[X„_i,Xi] n = 3,4,5,... (22) 

n — 2 

Y^:= ^ -[Y^_,,X^] m = l-l/z,2~l/z,3-l/z,... (23) 

m — 1/2 — 1 

As shown in Appendix A, any set of generators, which is contructed in such a way, 
satisfies all commutation relations in Eqs. ([1]) and ([2]) automatically. Therefore, we 
can conclude that any representation is fully determined by a single linear operator, 
namely, A2- 



5. Representation by integral kernels 

Given that the whole representation is determined by the linear operator A2, one has 
to specify the most general form of this operator under the constraints flT6|) . To this 
end it is convenient to represent An acting on some function in form of a convolution 
integral 

[A0](t,r) = jdt' jdr'An{t,r,t',r')<P{t',r') (24) 

with a kernel An{t,r,t',r') which can be thought of as the 'matrix elements' of An- 
Likewise, the operators Bm, which appear in the generator Ym, can be written as 
convolution integrals 

[i3^0](t,r) = jdt' jdr'Bn,{t,r,t',r')(j){t',r') (25) 

with a kernel Bm{t,r,t',r'). The local contributions of these kernels, which appear as 
ordinary differential operators in the LSI representation, correspond to Dirac (5-functions 
and their derivatives. However, in general the kernel may be non-local in space and 
time, including the recently discovered non-local representations involving fractional 
derivatives as special cases. 
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Because of the commutation relation [X„,Xo] = nX^ the kernel itself is a 
generalized homogeneous function under anisotropic dilatation by a factor 6 > 0: 

Anit/h, r/h^'\ t'/b, r'/b^/') = b^/'+^-''A^{t,r,t',r') . (26) 

Moreover, the kernel A2 is constrained by [[X2,X_i],X_i] = in Eq. flT6l) . tantamount 
to [[^2, dt], dt] =0, which implies that A2 has to be of the form 

^2(t,r,t',r') = ^(t + t')Kit-t',r,r')+L{t-t',r,r'), (27) 

where K and L are functions which depend only on the time difference t — t'. Because 
of Eq. (!26l) . they are generalized homogeneous functions as well: 

''inr'W'w) -'^"-'L(t-^,ry). (29) 

Hence any representation of the LSI algebra is determined by two homogeneous time- 
translation- invariant functions K and L. 



Because of Xi = |[X2,X_i] the kernel Ai = — |[^2,<9t] is essentially the temporal 
derivative of A2, i.e. 

Ai(t,r,t',r') = K{t-t',r,r'). (30) 

Because of Eq. (|2]), the operators Bm are given by 

Bk-i/z = ^^[Ak,dr], k = 0,1,2,... (31) 

so that the kernel of the shear generator Yi^i/z reads 

i?i_i/,(t, r, t', r') = id + |-) K{t - t', r, r') . (32) 

Again this kernel is a generalized homogeneous function, i.e., 

Bi-i/.(^, ^, p ^) = b-'^B,_y.{t,r,t',r'). (33) 

Once all kernels have been determined, one can easily check them by testing the 
commutation relations [X„,X_i] = (n -f l)X„_i and [y„,X_i] = (n + l/z)Yn-i, which 
can be translated into the differential equations 

(- + — ) An{t, r, t', r') = {n+ 1) A„_i(t, r, t', r') , (34) 

(- + — ) Bnit, r, t', r') = {n+ l/z)B„.^ {t, r, t', r') . (35) 
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6. Two-point correlation functions 

Let us now investigate the properties of a correlation function 

C{t, r, t', r') = (01 (t, r)02(t', r')) (t > t') (36) 

of two quasi-primary fields 0i and 02- For simplicity let us assume that both fields carry 
the same scaling dimension xi = a;2 = a; and that the correlator respects causality, i.e., 
it is nonzero only for t >t' . Because of translational invariance in space and time this 
function depends only on the differences of the coordinates. Moreover, invariance under 
global scale transformations, as expressed by the condition XqC = 0, implies the scaling 
form 

C{t, r, t', r') = e{t- t') {t - t')-'^/^ $ { JSjy/, ) , (37) 

where $ is a scaling function and 0(t — t') is the Heaviside step function which accounts 
for causality. 

In the following we first address the problem how $ can be determined for a given 
representation of LSI generators. Then we consider the inverse problem, i.e., for a given 
correlation function we ask for a suitable representation of the generators. 

6.1. Integro- differential equation for the scaling function 

LSI is based on the postulate that correlation functions of primary fields are invariant 
under the action of the generators, i.e., X„C = YmC = 0. These conditions lead to 
integro-differential equations which determine the form of the scaling function $. It is 
important to note that in the LSI theory the representation of the generators acting 
on 01 and 02 are generally different. More specifically, it was argued that a correlator 
vanishes unless the two representations are related in a specific way, giving rise to 
so-called Bargmann superselection rules. For example, in the standard Schrodinger 
representation the 'mass terms' occuring in the generators are known to have different 
signs. For this reason we will work with two different representations, denoting the 
integral kernels acting on 0i and 02 by the superscripts ^^^ and ^'^\ respectively. 

Let us first consider the generator Xi, which acts on the two-point function by 

0=[XiC](t,r, t',r') 

= - (t% + t"dt^ + —dr + ^dr^ + M + ^) c{t, r, t', r') (38) 

V z z z z J 

dr"K^^\t-t",r,r")C{t",r",t',r') 
dr"K^^\t'~t",r',r")Cit,r,t",r"). 
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Inserting the scaling form fl37p and using the homogeneity condition (1281) one obtains an 
integro-differential equation. As we have assumed the correlator to be translational 
invariant, this equation is generally over-determined unless the kernels K'^^'^-' obey 
specific constraints, referred to as Bargmann superselection rules. A straight-forward 
calculation shows that one obtains an autonomous integro-differential equation for the 
scaling function $ if and only if the two kernels have the form 

K^^\t-t\r,r')= hr + r'f f{t-t\r-r') (39) 

+ -(r + r') g{t -t',r - r') + h{t -t\r - r') 
K^^\t-t\r,r')= -hr + r'f f{t' -t,r' -r) (40) 

+ -(r + r') g{t' -t,r' -r) - h{t' -t,r' -r), 

where f,g,h are certain functions which depend on only two parameters. Because of 
Eq. ( l28l) these functions are homogeneous, hence they obey the scaling form. 

9ir,0=r-'/'G{^r-'/^), (41) 

With these kernels and the scaling form ( 1371) the integral equation (l38l) turns into an 
integro-differential equation for the scaling function $(0: 

<1.'K) + 2|^^d|^'d^^-^[|(l-^)-='V'F(|^^) + (42) 

As can be seen, the kernel function H drops out, meaning that it does not influence the 
form of two-point functions. Knowing the kernels F, G, this is the integro-differential 
equation which determines the scaling function $ of the two-point correlation function. 

Turning to the generator of shear transformations ^i-i/^, one can show that the 
invariance condition [Yi_i/^C](t, r, t', r') = leads exactly to the same integral equation 
and thus does not provide any new information. 

6.2. Inverse problem 

Let us now consider the inverse problem which can be formulated as follows: For a 
given two-point function characterized by x, z, and $ we would like to determine 
appropriate kernel functions K^^'"^^ such that $ is a solution of Eq. (l42l) . Since these 
kernel functions determine all LSI generators through recursion relations, solving the 
inverse problem would mean to construct a representation that renders exactly a given 
two-point function. In the following we outline how this problem may be solved. 
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For simplicity let us restrict to the case G = H = 0, i.e. we want to determine F 
for a given $. With the substitution ^ -^ l^~^^^i the integro-differential equation fH2l) 
turns into a convolution product 
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After integrating by parts in ^2 one can carry out the integration over i^. Finally another 
Fourier transformation of the entire equation yields 
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(46) 



This equation holds for all /c G M. As $ and F are symmetric, let us restrict to k > 0. 
Substituting A = /i^/^fc one obtains an inhomogeneous Volterra integral equation of the 
first kind for F{X) 



with the kernel 

^{k,X) 



dA^(A;,A)F(A) 
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Hence for given x, z, and $ the solution of the inverse problem amounts to perform the 
following non-trivial steps: 



(i) Compute the Fourier transform of $ in Eq. 

(ii) Plug $ into Eq. (HH]) . compute the kernel \l/(/i;, A), and solve the integral 
equation (H7|) to obtain F. 

(iii) Compute the inverse Fourier transform of F, plug it into Eqs. fHTl) and fl39|) and 

determine the kernel K. This kernel establishes the representation of Xi. 
(iv) Determine the LSI generators by recursion, setting e.g. L = 0. 
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7. Discussion 

The success of LSI applied to various exactly solvable system raised the hope that 
this theory might be a generic feature of anisotropic scale-free phenomena, including 
systems which cannot be solved exactly. To verify this expectation, various authors 
measured two-point autocorrelation and response functions numerically and compared 
them with the predictions of LSI using a suitable representation of its generators. For 
some systems, most notably the kinetic Ising model and the contact process, deviations 
were found, leading the authors to the conclusion that those systems are probably not 
invariant under local scale transformations. 

However, such a conclusion would be premature. The results of the preceding 
section lead to the conjecture that for any physically meaningful two-point function 
with arbitrary x and z it is possible to construct a representation of LSI generators 
which precisely reproduces the desired two-point function. On the one hand, this means 
that numerical discrepancies do not necessarily falsify LSI, instead they could also come 
from choosing the wrong representation. On the other hand, the space of possible 
representations is so huge that LSI itself can probably not predict the form of two-point 
functions, hence on this level it cannot be used to set up a classification scheme of 
anisotropic critical phenomena. 

For the sake of simplicity the present work was restricted to stationary situations 
in 1+1 dimensions, where correlation functions depend only on differences of the 
coordinates. It would be interesting to investigate recent applications of LSI to 
ageing [4,5,9,19,21,33-38] along similar lines. Here the algebra is relaxed by giving 
up time-translational invariance, probably reducing the predictive power of the theory 
even further. 

Why is LSI less predictive than conformal invariance in two dimension? In my 
opinion this issue can be traced back to different symmetry properties. Conformal 
invariance [2] generalizes global dilatations combined with rotations to a local symmetry. 
Similarly, LSI generalizes global dilatations combined with shear transformations to 
a local symmetry. However, rotations and shear transformations are very different 
in character. In the case of conformal invariance, an isotropic equilibrium model is 
expected to be rotationally invariant hy itself. Contrarily, as demonstrated in Fig. [21 
an anisotropic process is not automatically invariant under shear by itself, rather this 
invariance has to be established manually by adding suitable terms to the generators. At 
this point the theory of LSI requires an input of extra information which is not needed in 
conformally invariant systems. The message of this paper is that the ambiguity caused 
by this additional information reduces the predictive power of LSI. More specifically, it 
is suggested that this information can be expressed in terms of functions with a single 
parameter, providing so many degrees of freedom that any two-point function can be 
reproduced by the theory. 
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The present findings make it plausible why LSI was apphed successfully to many 
exactly solvable systems while it continued to fail for certain non-integrable systems. 
In exactly solvable systems one usually arrives at a partial differential equation which 
determines the two-point function. This allows one to derive suitable LSI generators 
in a closed form. For a non-integrable system such as directed percolation [26] the 
scaling function $ is a complex object which involves loop corrections to all orders of 
the underlying field theory. The corresponding LSI generators may exist, but apparently 
it is impossible to write them down in a closed form. In such cases the attempt to guess a 
suitable representation and to confirm it numerically by comparing two-point functions 
is likely to fail. 

The results of the present work do not rule out that LSI might have some predictive 
power on the level of three-point functions. In my opinion this is one of the key issues 
to be addressed in the future. 
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Appendix A. Consistency of the contruction scheme 

In this appendix it is shown that the generator X2, constrained by Eq. (fT6!l . determines 
aU generators iteratively in such a way that all commutation relations are satisfied. 

A.l: Commutators [X„,Xm] = {n — m)Xn+m: 

To prove these commutation relations by induction, we first notice that 

[Xo, X_i] = X_i, [X2, X_i] = 3X1, (A.l) 

[X2, Xq\ = 2X2, [Xi, X_i] = 2X0 . 

The first relation is always fulfilled, while the second one was used to define Xi in 
Eq. (1211) . The third and the fourth relation have been used as constraints for X2 in 
Eq. ( TT6l) and thus they are satisfied as well. Moreover, the generators X^ have been 
constructed according to Eq. fl22|) . hence the commutators 

[X„,Xi] = (n-l)X„+i n = 2,3,... (A.2) 

are valid by construction. Anchored at these relations, the remaining commutation 
relations [X„,Xm] = {n — m)Xn+rn can be proven by induction. To this end let us 
assume that these commutation relations hold for n + m = N — 1 and show that they 
also hold for n + m = N: 

[An,,Xm] = -[[X„_i,Xi], Xm] 

= ^^([[X„,Xi],X„_i] + [[X„_i,X„],Xi]) (A.3) 

apply 1I22I 1 induction 

= ^( ("^- l)[Xm+i,X„_i] + {n-m-l) [X„+m_i,Xi]) 

n — Z \ -^ V '^ 

apply ^ 

= —^{i^- l)[X^+i,X„_i] + {n-m- l)(r2 + m-2)X„+^]j 

This is again a recursion relation with fixed N for an inductive step from (n — 1, m + 1) to 
{n,m). More specifically, if the commutation relation [Xm+i,X„_„i] = (m — n + 2)Xm+n 
is known to be valid, the above equation implies that the commutator 

{1 — m){n — m — 2) — {m — n + l){m + n — 2) 
n — 2 
= {n- m)Xn+m (A.4) 

The same recursion works also in different direction as an inductive step from (n+1, m — 
1) to {n,m). This twofold recursion scheme allows one to check all commutators of the 
form [X„,Xm] = (n — m)Xn+m iteratively. 
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A2: Commutators [X„, y^] = {n/ z — m)Yn^rn'- 

To prove these commutations relations, we first show by induction that 

n -\- 1 

[X„,F_i/,] = ^F„_i/,. (A.5) 

The induction is anchored at ra = — 1 since [X.ijYli/^] = [—dt,—dr] = 0. Assuming 
that the commutation relations are satisfied for n — 1, i.e. 

Tl 

[A„_i,r_i/,] = -F„_i_i/, (A.6) 

the same relations holds for n because of 

[Xn,Y_y,] = -1^ [[A„_i,Xi],r_i/^] (A.7) 



n 



i^([Xi, [F_i/„X„_i] +[X„_i, [Ai,y_i/,]] 





1 ( 


apply (|A.6|| 

[-^l> ln-l-l/z\ H 

apply 1(23 

J z\z J 


apply II23II 

Xn-1, Yl-l/z] 
apply (IA.6II 


n 


-2\ 

z\z 



n — 2 z 

Next, let us assume that the commuation relations 

[X„, Fm_i] = {njz - m - l)F„+m-i (A. 8) 

hold for all n and a given m, anchored at m = — 1/2; by Eq. flA.Sp . Then we can prove 
the remaining relations by a another induction: 

[X„, Ym\ = -^ [X„, [Ai, y^_i]] (A.9) 

- + 1 — m 



y— — [ [ J^m- l i^nL ^l] + [[A„, Ai], l^_i] j 

apply IIA.8I 1 apply JlJ 

1 / Ti 

_(_ _ ^ + 1) [F„+„_i, Xi] +{n - 1) [A„+i, y^_i] 

- + 1 — m V z ^ ..^ / V ^^ /, 

apply ^^ apply IIA.8I I 

[f - m + 1 Vi - m - n + 1) + (n - 1) [^ - m + 1 



V 

1 , T -" n+m 

- + 1 — m 



m] Y„ 
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Appendix B. The most important representations 

This appendix demonstrates that the most important representations of LSI derived so 
far can be described within the unified framework of generating kernel functions. 

a) Schrodinger representation 

The standard Schrodinger for diffusive systems with z = 2 is defined by [39-41] 

Xn= - t-^% - ^t-rdr - |(n + ir - ^^^^^A^r-V^ , (B.l) 

Y^= - r+^/'S, - (m + l/2)Mt"'-'/\- , (B.2) 

where Ai is the so-called mass parameter. As can be verified easily, this representation 
is generated by the kernel functions 

K{t-t',r,r') = ^Mr^6{r-r')6{t-t'), L{t - t' ,r,r') = . (B.3) 

b) Local representation type (i) for arbitrary z 

This representation is given by [1] 

X„ = - r+i^i - '^^t'^rdr - -(n + l)r - I^il^dlHsior-V^ ,(B.4) 

z z 2 

n-1/. = - t'^dr - ^kB,ot'-'r'-' (B.5) 

and reduces to the Schrodinger representation for z = 2. The corresponding kernel reads 

K{t-t',r,r')=Bior'6{r-r')6{t-t'), L{t - t',r,r') = . (B.6) 



b) Extended local representation type (ii) for z = 2 

The representation [1] 

A„ = - t-^% - -^t-rdr - -(n + l)r (B.7) 

^(^ + 1) ^ -n-i 2 {n^~l)n 2 4 

;; -Diot r ±>20t r 

2 o 

n-i/. = - f'dr - 2kBwt'-'r - h{k - l)B2ot'-\^ (B.8) 

involves both kernels K and L: 

Kit - t', r, r') = Bi^Sir - r')5{t - t') , (B.9) 

L{t-t',r,r) =B2or^S{r-r)6{t-t'). (B.IO) 
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c) Local representation type (iii) for z = 1 

The 'conformally invariant' representation [1] 

X, = - r+i^i - A^,' [{t + Aior)"+i - r+i]9, 
- (n + l)a;t" - 



^ + '^^°[(t + A,orr-r] 



A 



10 



k 



Yk-i = - (t + Aior)'^9, - -5io(t + Aior 



»fc-i. 



corresponds to 



Kit - t', r, r') = 5{t - t') Ai^d'^r - r') + Bi^rdij- - r') 



L{t - t', r, r') = 6{t - t') Al^r'^d'ij- - r') + -Bi^AiySij- - r') 
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(B.ll) 

(B.12) 

(B.13) 
(B.14) 



d) Temporally nonlocal representation for arbitrary z 

Tliis representation, called type I in Ref. [30], is given by 



Xi = -t'dt--trdr 1 (B.15) 

z z 

- (/3 + -fyd^,-' - 27(2 - z)rdl-' - 7(2 -z){l- z)d;' 

Yx-Mz = -tdr-{(3 + ^)zrdl~' - 7^(2 - ^)9^^ . (B.16) 

It is non-local in space and involves fractional derivatives, depending on z even with 
negative powers. Fractional derivatives can always be expressed by integral kernels. As 
one can see, this would determine the kernel K while the kernel L vanishes. 



e) Spatially nonlocal representation for arbitrary z 

Another recently discovered representation, called type 11 in Ref. [30], is given by 

Xi = - t'^dt - -trdr -—t- ar'^df'-^ (B.17) 

z z 

Yx-xi-, = -tdr- azrdl'"'^ . (B.18) 

As it involves temporal fractional derivatives, this representation is non-local in time. 
Again this would determine the form of the kernel K while the kernel L vanishes. 
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